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ABSTRACT 

We investigate the star formation and dust extinction properties of very luminous infrared 
galaxies whose spectra display a strong US line in absorption and a moderate [On] emission 
(e[a] spectrum). This spectral combination has been suggested to be a useful method to 
identify dusty starburst galaxies at any redshift on the basis of optical data alone. We compare 
the average e(a) optical spectrum with synthetic spectra that include both the stellar and the 
nebular contribution, allowing dust extinction to affect differentially the stellar populations 
of different ages. We find that reproducing the e(a) spectrum requires the youngest stellar 
generations to be significantly more extinguished by dust than older stellar populations, and 
implies a strong ongoing star formation activity at a level higher than in quiescent spirals. A 
model fitting the optical spectrum does not necessarily produce the observed FIR luminosity 
and this can be explained by the existence of stellar populations which are practically obscured 
at optical wavelengths. Models in which dust and stars are uniformly mixed yield a reddening 
of the emerging emission lines which is too low compared to observations: additional foreground 
reddening is required. 
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1. Introduction 



The effects of dust extinction are crucial for interpreting the spectra of nearby and distant galaxies, 
particularly after that new facilities from the ground (e.g. SCUBA on JCMT) and from space {ISO) 
have established that the IR emissivity of galaxies was much enhanced in the past with regard to what is 
typically observed in local galaxy populations (Elbaz et al. 1999; Barger et al. 2000, Franceschini et al. 
2000, Ivison et al. 2000, Small et al. 2000). Further evidence about the very important role played by dust 
reprocessing during the evolution of galaxies was provided by the detection of a diffuse far-IR background 
likely originating from distant and primeval galaxies, whose dominance over the optical background testifies 
hidden star formation at unexpectedly high rates (Puget et al. 1996; Dwek et al. 1998). 
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Waiting for future spectroscopic capabilities in the far-IR which will be offered by a variety of missions 
(SIRTF, FIRST, and NGST), only indirect clues about the relation between dust and stellar population 
emissions can be gained from optical spectra. 

In the local universe, a large fraction of the FIR-luminous galaxies are characterized by a peculiar 
combination of spectral features in the optical: a strong KS line in absorption (EW> 4 A) and a moderate 
[On] emission (EW< 40 A) (Poggianti & Wu 2000, hereafter PWOO). Galaxies with this type of spectra 
were named "e(a)" galaxies f| and were found to be quite numerous in the cluster and field environments at 
z = 0A- 0.5 (Dressier et al. 1999, Poggianti et al. 1999, hereafter P99). The equivalent width of their US 
line exceeds that of typical, quiescent spirals at low-z and their low [On] /Ha ratios are consistent with the 
emission line fluxes being highly extincted by dust (P99, PWOO). 

Interestingly, the combination of moderate emission lines and strong early Balmer lines in absorption 
has been subsequently found also in the first spectroscopic surveys of the high-redshift galaxy populations 
detected by ISO (Aussel 1999, Flores et al. 1999). The association between the e(a) signature and the 
FIR-luminous galaxies seems to be confirmed at any redshift and suggests that the e(a) spectra - despite 
their moderate emission lines - actually belong to highly obscured, starburst galaxies. One of the most 
important pieces of information that is still missing is the physical origin of this type of spectra: what star 
formation and dust properties produce this peculiar spectral combination? P99 and PWOO proposed that 
the e(a) spectrum originates from "selective extinction" in a galaxy where the youngest stars in the HII 
regions are much more dust-embedded than the older stellar populations, which are more widely distributed 
throughout the galaxy. Such an age-dependent obscuration would cause the emission lines to be highly 
extincted (hence producing the low blue-to-red emission line ratios observed) and would allow the older, 
A-type stars to dominate the integrated spectrum at A ~ 4000 A (thus explaining the strong US line in 
absorption) . 

The existence of a selective extinction in galaxies is well known, but it is usually ignored when 
dealing with the global effects of dust on the spectral energy distribution of galaxies. Stars are expected 
to spend the very beginning of their evolution deeply embedded in dusty environments, later drifting 
away from or dispersing the molecular clouds where they were born (Calzetti et al. 1994 and references 
therein, Granato et al. 2000). As discussed in PWOO, an age-dependent dust obscuration could explain the 
numerous observational studies which measure discrepant extinction values when using different spectral 
ranges/features (Israel & Kcnnicutt 1980, Viallefond et al. 1982, Caplan & Deharveng 1986, Fanelli et al. 
1988, Bohlin ct al. 1990, Keel 1993, Calzetti et al. 1994, Veilleux et al. 1995, Lancon et al. 1996, Mayya & 
Prabhu 1996, Calzetti 1997,1998, Mas-Hesse & Kunth 1999). Moreover, Zaritsky (1999) found a clear trend 
of decreasing mean extinction with declining stellar effective temperature (i.e. stellar age) in the Large 
Magellanic Cloud. 

So far, no attempt has been made to model the whole e(a) spectrum in a quantitative way. The aim 
of this work is to test whether the average e(a) spectrum can be modelled as a current starburst assuming 
a selective dust extinction that increases as the age of the stellar population decreases. In the following we 
employ a spectrophotometric model that includes both the stellar component and the thermal emission of 



^We stress that these are not the so-called "E-l-A" (or k+a) spectra which by definition have no detectable emission 
lines (Dressier & Gunn 1983), at least at the level of the [On] detection limit of the spectroscopic surveys of distant 
galaxies (typically 3 A). Nevertheless, it is possible that some of the k-l-a spectra are e(a) galaxies with an extremely 
low [On] (P99, Smail et al. 1999). 
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the gas ionized by young stars. Such a model wiU be considered successful only if it is able to reproduce 
simultaneously the shape of the optical observed continuum (3600-6900 A) and the strength of the most 
important emission/ absorption lines: [Oii]A3727, US at A = 4101 A, H/3 at 4862 A and Ha at 6563 A. 

The plan of the paper is as follows: in §2 we present the spectroscopic catalog and the average 
observed e(a) spectrum and in §3 we describe the ingredients and the structure of our model. The results 
are discussed in §4, where we examine starburst models with selective extinction (§4.1), post-starburst 
models (§4.2) and models with constant extinction (§4.3). Our main conclusions are summarized in §5. 



2. Observations 

The observed e(a) spectra are taken from the complete spectroscopic catalog of Very Luminous Infrared 
Galaxies {\og[Lipi/ Lq) > 11.5, VLIRGs) presented in Wu et al. (1998a, b). Among these, we selected those 
19 galaxies with a secure e(a) spectrum and a spectral coverage extending below the [On] line and beyond 
the Ha line. These are nuclear spectra covering at least the central 2 kpc of the galaxy (see also §4). None 
of these galaxies is an AGN (either Seyfertl or Seyfert2) according to the emission line ratios (PWOO). The 
average e(a) spectrum is shown in Fig. 1 and the equivalent widths of the main lines are listed in Table 3. ^ 

As mentioned in §1, moderate emission lines and strong early Balmer lines in absorption have been 
found also in the spectra of the high-redshift IR populations. In fact, the average spectrum of the distant 
starburst galaxies detected at 15 /im by ISO in one of the fields of the Canada- France Redshift Survey 
(Flores et al. 1999) turns out to be very similar to the average e(a) spectrum. The comparison is presented 
in Fig. 1 and shows that ~ within the wavelength range in common (-3700-5100 A) - both the continuum 
shape and the strength of the [On], H(5 and H/3 lines are alike. 
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Fig. 1. A rest-frame comparison of the average e(a) spectrum of local Very Luminous Infrared galaxies (thick line) and 
the spectrum of the ISO starburst galaxies at < ^ >~ 0.5 (thin line). The latter is taken from Fig. 12 in Flores et al. (1999) 
and is the average spectrum of 5 ISO-detected galaxies at z < 0.7 whose spectral energy distributions at visible, near-IR, MIR 



^Spectra have not been corrected for foreground Galactic extinction, which is negligible (mean Galactic E(B-V) 
= 0.04) in comparison with the intrinsic extinction of these galaxies (Wu et al. 1998b). 
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and radio wavelengths resemble those of highly reddened starbursts in the local universe (Schmitt et al. 1997); it has been 
normalized to the VLIRGs spectrum over the wavelength range in common. The spectral resolution is 10 A (Wu) and 40 A 
(Flores). 



3. Models 



Instead of employing a highly complex tool (i.e. a complete spectrophotometric model), we preferred 
to develop and use a "simplified model" . This model retains all the basic physical ingredients required 
by the complexity of an integrated spectrum, but includes only those stellar populations that are a priori 
essential because they are known to affect the spectral features that we wish to reproduce. 

The integrated model spectrum has been generated as a combination of 10 stellar populations of 
different ages, listed in Table 1. Each stellar generation is born with a Salpeter IMF between 0.1 and 
100 AIq. The ages of the 10 populations have been chosen considering the evolutionary timescales 
associated with the observational constraints: the youngest generations (10^,3 • 10^,8 • 10^,10^ yr) 
are responsible for the ionizing photons that produce the emission lines; the intermediate populations 
(5 • 10*", 10^, 3 • 10^, 5 • 10®, 10^ yr) are those with the strongest Balmer lines in absorption^, while older 
generations of stars (named "old" in Table 1 ) have been modelled as a constant star formation rate (SFR) 
between 1 and 12 Gyr before the moment of the observation and can give a significant contribution to the 
spectral continuum, hence affecting also the equivalent widths of the lines. This approach allows us to a) 
reduce the number of parameters to a minimum: a complete spectrophotometric model finds the integrated 
spectrum summing over typically 50 single stellar populations of various ages, instead of 10; b) explore 
also bursting, discontinuous star formation histories which more likely resemble the irregular histories of 
"real" starburst galaxies. Adopting as star formation history an analytic function + a single or multiple 
burst would force us to restrict the parameter space investigated according to some a priori assumptions; c) 
easily implement different extinctions for stellar populations of different ages, interpreting the results in a 
straightforward manner. 

The composite (stars+gas) spectrum of each single generation has been produced with the 
spectrophotometric code of Barbaro & Poggianti (1997, BP97), that includes both the stellar component 
and the thermal emission of the ionized gas. This model is fully described in BP97; here we only summarize 
its salient characteristics. Between 3500 and 7500 A, the stellar model is based on the stellar spectral 
library of Jacoby et al. (1984) which has enough resolution to study absorption features, such as the stellar 
lines of the Balmer series in absorption. The ionizing flux of the young stellar populations is computed 
using stellar atmosphere models (Kurucz 1993); from this, the luminosities of the nebular hydrogen lines 
(Balmer series) are derived assuming case B recombination (Osterbrock 1989) while the strength of other 
lines - such as [Oii]A3727 ~ is calculated using HII region models (Stasinska 1990). When non- negligible 
(for ages < a few 10'' yr), the nebular contribution (lines + nebular continuum) is added to the stellar 
component to give the total spectrum. For this project, both the stellar and the nebular components were 
found assuming solar metallicity. The model spectra were broadened with a Gaussian degrading the original 
model resolution (~ 4 A) in such a way to match the resolution of the observed spectra (~ 10 A). 



^The EW(H(5) is always > 4 A in the spectra of these intermediate populations. 
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Each single generation (SG) is assumed to be extincted by dust in a uniform screen according to the 
standard extinction law of the diffuse medium in our Galaxy (Ry = Av/E{B — V) = 3.1, Cardelli et al. 
1989) unless otherwise stated (see §4). The extinction value E(B-V) is allowed to vary from one stellar 
population to another and the extincted spectral energy distributions of all the single generations are added 
up to give the total integrated spectrum. 

Within a chosen star formation scenario (see §4), the best-fit model was obtained by minimizing the 
differences between selected features in the observed and model spectra. A merit function was constructed 
considering TV = 12 features: the equivalent width of four lines ([OII]A3727,H(5, H/3 and Ha) and the relative 
intensities of the continuum flux in eig ht almost featureless spectral windows (3770-3900A, 4020-4070A, 
4150-4250A, 4600-4800A, 5060-5150A, 5400-5850A, 5950-6250A and 6370-6460A): 



where the quantities M^, O;, E^, Eq and Wi refer to the value predicted by the model, the observed 
value, the accuracy, the minimum error and the weight assigned to the i feature, respectively. An accuracy 
of 1% ad 10 % was adopted for the flux in the continuum bands and for the equivalent widths, respectively. 
These are comparable to the observational uncertainties if we consider that the error on the relative flux 
calibration amounts to a few percent, the typical error in the measured equivalent widths is between 10% 
and 20% and we are dealing with an average spectrum.^ 

All quantities were given a unit weight but for the 3770-3900 A continuum and for the [OIIJA3727 line, 
to which a lower weight was assigned (0.4 and 0.5, respectively). This choice is due to the enhanced flux 
calibration uncertainty of Jacoby's spectra below 4000 A and to the fact that, unlike the Balmer emission 
lines, the intensity of the [OIIJA3727 line is not related in a simple way to the star formation rate. In fact, 
the [On] line is sensitive to several other factors such as the hardness of the ionizing spectrum and the local 
conditions in the nebula, the density and the metallicity of the gas. We find that fitting the observed [On] 
strength requires a line intensity higher than the one predicted by models of solar metallicity ^ unless the 
presence of hot Wolf Rayet stars is included. At ages between 3 Myr and 5 Myr these stars may enhance 
the ratio [OII]A3727/Ha by a factor of about three (Garcia Vargas et al. 1995), hence we chose to take this 
factor into account increasing the [On] flux of the 3 Myr old single generation. 

The merit function was minimized with the method of the simulated annealing (routine AMEBSA in 
Numerical Recipes), which is particularly suited to avoid the rapid convergence to a local minimum. The 
free parameters (the relative intensity of the SGs and the extinction values) were let to vary within specified 
ranges defined by the kind of evolutionary scenario considered. 

In principle, the observed FIR/optical luminosity ratio could be treated as a further constraint of the 
model, R hence, in order to predict the FIR emission, we computed a new set of single stellar populations 



^In addition, we checked that assuming a lower accuracy (4% and 20% in the continuum and in the equivalent 
widths, respectively) does not modify our main conclusions. 

^We note that a metallicity below solar produces a stronger [On] line than the solar case, but causes also an 
enhancement of the [OIII]A5007 line which instead is quite weak in the e(a) spectrum. 

^The FIR luminosities given by W98a are computed from the IRAS 60 ^m and 100 ^m fluxes as Fir — 



N 



Ei + Eq 




(1) 



i=l 
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by extending the Jacoby stellar spectra below 3510 A and above 7400 A using Kurucz (1993) stellar 
atmosphere models, with a procedure similar to that adopted by Longhetti et al. (1999). Once the 
parameters of a best-fit model were determined, the FIR luminosity was evaluated as the difference between 
the total unextinguished spectrum and the final (extinguished) spectrum. 



4. Results 

In this section we show the results of some selected models whose best-fit parameters are listed in 
Tables 1 and 2. For each model, the left column (SF) gives the star formation rate (mass per year, after 
having normalized to SFR=1 yr^^ in the old population) and the right column gives the color excess 
E(B-V). The visual comparison of these models with the average e(a) spectrum is shown in Figs. 2, 3 and 
4 while Table 3 presents the equivalent widths of the four lines considered, the FIR to optical luminosity 
ratio and the current star formation rate of each model. 

Different evolutionary scenarios have been examined in order to address the following questions: 

a) Can the star burst, selective extinction scenario explain the e(a) spectra? 

b) What constraints can be placed on the burst duration and strength? 

c) Can the post-starburst scenario explain the e(a) spectra? 

d) Can the e(a) spectrum be fitted by models with a constant value of extinction? 



4.1. Can the starburst selective-extinction scenario explain the e(a) spectra? 

At first we consider models in which the star formation rate is roughly constant (within a factor ~ 1.5) 
during the current starburst and much lower before that (models A and B in Table 1 and Fig. 2). The 
extinction varies with the age of the population, being significantly higher at younger ages. 

Model A is an example of a starburst that began 2 • 10^ yr ago. The quality of the fit can be assessed 
from the difference between the model and the observed spectrum in Fig. 2: both the line strengths and 
the highly reddened continuum are well reproduced, except for the continuum level around 5000 A where 
anyway the discrepancy is only at the 4 % level. This problem is common to all the models, regardless of 
the chosen set of parameters (Figs. 2, 3 and 4). 



1.75[2.55Seo + l-OlySioollO ^'^ Wm ^ and are approximately total FIR luminosities in the range 1-1000 fxm. 
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A(A) log(AGE) (yr) 



Fig. 2. Left: Comparison between tlie average observed spectrum of the e(a) galaxies (thick hne) and the spectrum of the 
starburst models listed in Table 1 (thin line) normalized at 5500 A. The emission lines of the models have been assigned a 
Gaussian shape of (arbitrarily) small width, hence their equivalent width, not their shape, is relevant for the fit (see Table 3). 
The difference between the model and the observed spectrum is also shown. When plotting this difference we have not 
considered the emission lines which we did not attempt to model ([OIII]4959 and 5007 and [Ol]6300) and the Nal D line at 
5893 A which is mostly of interstellar origin. The vertical segment on the right side represents the 1 a error in Ha. The 
difference in the flux of the lines has been plotted as a Gaussian with a width equal to that of the observed line. Right: SFR 
(histogram, normalized to =1 in the old population) and E(B-V) (dots) of the models whose spectrum is shown in the left 
panel. The observed FIR/V ratio is =88.0 (Table 3). 



Model A can only account for about 1/3 of the observed average FIR luminosity {Lfir/Lv — 28.9 
instead of 88.0). There can be several reasons for this disagreement: first of all, it should be kept in 
mind that the observed spectra refer only to the central region of the galaxy as it appears in the opticalF 



^In a typical FIR galaxy, about 2/3 of the observed Hq luminosity come from regions more than 1 kpc from the 
nucleus (Armus et al. 1990). Neverthless, considering that the color excess measured from the Balmer decrement 
strongly decreases with radius (Veilleux et al. 1995, Kim et al. 1998) and that the primary mechanism for ionizing 
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However, also the IR emissivity is usually concentrated in the central regions of luminous infrared galaxies 
(see e.g. Soifer et al. 2000, Kennicutt 1998 and references therein), therefore aperture effects are not likely 
to be responsible for the remaining 2/3 of FIR emission. 

Moreover, there might be star forming regions that are completely obscured at optical wavelengths. 
This is known to be a common situation in Luminous Infrared galaxies and a spectacular example is the 
image of the Antennae galaxies (Mirabel et al. 1998). The most intense starburst in this merging system 
takes place in a region that is inconspicuous at optical wavelengths and coincides neither with any optically 
bright region nor with the dark lanes observed in the optical image. Such a starburst would give no 
contribution to our modelled spectrum, while accounting for a significant fraction of the FIR luminosity; 
therefore, our estimate of the current star formation rate relative to the old stellar generations would be 
correspondingly underestimated. 

Model B is an example that describes a starburst with a high FIR/V ratio due to some practically 
obscured regions. As in model A, the starburst began 2 ■ 10* yr ago, but both the SFR during the burst 
and the extinction of the two youngest stellar generations are higher than in A. The fit of the observed 
spectrum is again remarkable, with a small difference in Ha being well within the 1 a error (Fig. 2). Thus, 
model B can simultaneously account for the optical properties and the FIR luminosity. 

For a standard Salpeter IMF in the range 0.1 — 100 Mq, in model B about 60% of the total mass in 
stars is formed during the burst. Assuming a top-heavy IMF during this phase would substantially reduce 
the mass fraction in young stars: for example, with a Salpeter IMF lower mass limit = 1 Mq, the starburst 
in model B forms about 40% of the total stellar mass. ^ 

Within the starburst scenario with selective extinction, we find an excellent agreement between the 
model and the observed spectrum as far as both the continuum shape and the strength of the lines are 
concerned. A model that fits the optical spectrum doesn't necessarily reproduce the observed FIR/optical 
ratio but the modelled and observed ratios can be reconciled taking into account young stellar populations 
mostly obscured by dust at optical wavelengths. 



4.1.1. Non-standard assumptions regarding the dust: extinction law and relative distribution 

dust-stars 



the emission-line gas outside of the nucleus could be shock-excitation from outflowing winds instead of heating by 
young stars (Armus et al. 1989), the majority of the intrinsic Ha luminosity (after correcting for extinction) related 
to star formation should arise from the central regions. 

*The spectral energy distribution of a young (e.g. < 2 x 10** yr) stellar generation with an IMF between 0.1 and 
100 Mq is identical to that of a population with an IMF truncated at low masses (e.g. rriinf = IMq) because the 
contribution of stars with M < IMq is negligible at these ages. Hence, a model with a standard IMF (any of those 
presented in this paper) gives the same results of a model with a truncated IMF during the burst with the same set 
of E(B-V) values and the same star formation rate in stars more massive than 1 Mq: the only difference between 
these two models will be the total star formation rate during the starburst, with the truncated model missing stars 
below 1 Mq. 
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The models described so far have been computed assuming a standard Galactic extinction law and a 
dust screen around each stellar generation. It is interesting to investigate whether a satisfactory fit can be 
found imder other assumptions regarding the dust properties. 

In model C we consider a 2 • 10* starburst and an extinction law with Ry = Av/E{B — y) = 5 as 
observed towards some dense clouds in our Galaxy (Mathis 1990). The fit to the e(a) spectrum is still 
acceptable but the fit of the continuum around the Ha line worsens and there is a small difference in the 

Ha emission at the ^ la level (Fig. 2). In this case the observed FIR/optical ratio is fully reproduced even 
with E(B-V) values generally lower than in B because, for a given E(B-V), the optical flux is much more 
extinguished. 

Finally, we modified the assumptions regarding the relative distribution of dust and stars: instead of 
adopting a dust screen, we examined the case of internal dust, in which dust grains and stellar populations 
are uniformly mixed, for a standard extinction law {Ry = 3.1). In this case, the relation between the 
observed intensity and the intrinsic intensity 7° is the following: = Ix^~T~(^J~' ^^i^re rsc(A) is 
computed according to Calzetti et al. (1994). Interestingly, it was impossible to find a good fit to the 
observed spectrum. This is due to the fact that in the internal dust models increasing the obscuration (i.e. 
the intrinsic optical depth) does not yield a corresponding increase in the reddening of the spcctnmi: the 
latter saturates to an asymptote, and the asymptotic reddening value (E(B-V)~ 0.18)) is too low to be able 
to account for the observed emission line ratios (E(B-V)~ 1.1 from the observed Balmer decrement, W98a). 
This saturation of the mixed case was noted by Calzetti et al. (1994, Sec. 5. 3) and our asymptotic value of 
E(B-V) is in agreement with their asymptotic value of the measured Balmer optical depth (r^ = 0.24). This 
result can be easily understood considering that the most heavily reddened stars are those most heavily 
obscured, hence they cannot provide the dominant contribution to the emerging light, which instead must 
come from the outermost layers (Witt et al. 1992, Calzetti 1998): in the internal dust model this effect 
gives rise to a natural limit to the reddening. While the existence of regions where dust and stars are 
mixed can remove the contribution to the optical spectrum of entire stellar populations (those most deeply 
embedded), these regions do not appear to determine the characteristics of the emerging spectrum, which 
requires a higher reddening at least partly due to foreground dust. 



4-1-2. Burst duration and strength 

How long (short) can the burst be? It is interesting to investigate whether an e(a) spectrum (in particular, 
its US strength) requires a quite long starburst involving intermediate-age populations (10* — lO'' yr). 
Contrary to the most intuitive expectations, this is not the case. As an example of this, we present the 
good fit obtained for a starburst that began only lO'' yr ago (model D). In this case the strong H^ line 
is not produced by the stars born in the starburst event, but by the previous stellar generations. We find 
that starbursts as short as a few times 10^ yr and as long as 10^ yr are able to produce the e(a) spectrum 
as long as: a) the youngest populations (< a few lO'' yr) must be highly extincted and b) there must be a 
contribution of the "old" population, but this should not overwhelm the intermediate-age contribution at 
~ 4000 A. 

Model D has a FIR/V=62.8 therefore accounts for more than 2/3 of the observed FIR luminosity, but 
as for longer starbursts (models A and B), a good fit to the optical spectrum can be obtained with short 
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starburst models spanning a wide range in FIR/V ratio, simply modifying the amount of "hidden star 
formation" . 

As an extreme example of what we have just shown and to further illustrate the degeneracy in the 
physical solutions allowed by optical data, we consider a galaxy with a practically constant SFR during its 
whole evolution (model E) and we show that this can still fit the e(a) spectrum reasonably well - with a 
slight mismatch in Ha and in the continuum around it (Fig. 2) - if we adopt high extinction values for the 
young populations and a quite strong extinction of the old population (E(B-V)=0.5). It is worth stressing 
that the E model does not represent a "spiral-like" star formation history, because in all but the latest 
spiral types the SFR at recent times is known to be lower than the mean SFR during the previous evolution 
(Kennicutt 1998). We also tried to fit the e(a) spectrum assuming a spiral-like history in which the SFR 
of the old population is at least 1.5 times the SFR of the young and intermediate populations. This model 
cannot fit the observed spectrum unless we allow the extinction of the "old" population to be unphysically 
high, higher than that of the intermediate-age populations and much higher than the average observed even 
in the HII regions of nearby spiral galaxies (about 1 mag at Ha, corresponding to E(B-V)'~ 0.4). 

Although model E is able to reproduce the optical spectrum, its FIR/V is a factor of 10 lower than 
what is observed, and in this case (i.e. for an approximately constant star formation) no young obscured 
population can be invoked to explain the missing FIR flux; any additional current star formation would 
qualify such a model as a star burst. 
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POST-STARBURST MODELS 




4000 5000 6000 7000 6 7 8 9 10 

A(A) log{AGE) (yr) 



Fig. 3. Results of the post-starburst models listed in Table 2. Symbols as in Fig. 2. In the lower panel, the 1 a errorbar 
in H/3 is shown on the left side. 



4.2. The post-starburst hypothesis 

Here we consider another family of models, in which the galaxy is seen after a strong starburst during 
the post-starburst phase (Table 2 and Fig. 3). As discussed in P99, the two spectral features that define 
an e(a) spectrum (a strong H(5 line and a moderate [On]) can be reproduced by dust-free models of a 
post-starburst galaxy with some residual star formation, but on the basis of the FIR properties of the e(a) 
galaxies and of their unusually low [On] /Ha ratios, the same authors concluded that e(a) spectra more likely 
belong to starburst galaxies that are highly extincted by dustj^ Here we contemplate post-starburst models 
with a small residual star formation activity and including extinction by dust. This family of models shows 
a tendency either to underestimate the Ha strength or to overestimate H/3 (see models PSBl and PSB2 
in Fig. 3 and in Table 3): the emission in at least one of these two lines is more than 1 a away from the 
observed value. The reason for this is that the observed Ha/H/? ratio requires a high extinction of the lines 
(the extinction measured from this line ratio is E(B-V)~ 1-1), but such a high extinction in a post-starburst 
galaxy would produce too weak emission lines i.e. there isn't enough current star formation to give the 
observed strength. Although the discrepancy between the model and the observed H/3 strength in model 
PSB2 is not dramatic on an absolute scale (still less than 3 A), the post-starburst interpretation is not 
favoured even on the basis of the optical data alone (without taking into account the strong FIR emission) 
i/ the whole spectrum including the Ha line is examined. In addition, the FIR luminosity produced by these 



® The post-starburst interpretation has been mentioned by Flores et al. (1999) (see §2), who consider their data 
to be consistent with the last phases of a burst following strong star formation episodes, where the MIR emission is 
still high due to dust heating by intermediate-mass stars (A/ — 1 — 3AIq). 
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models is just a small fraction of the observed value and, similarly to model E, within the post-starburst 
scenario no young obscured population can be invoked to explain this discrepancy. 



4.3. A reference case: models with constant extinction 

So far we have only considered models in which the E(B-V) is allowed to vary from one stellar 
population to another. In this section we contemplate more "canonical" models which have a common value 
of extinction for all the stellar generations, showing that these fail to reproduce the e(a) spectrum. 

Having imposed the extinction to be the same for all the 10 populations and to be within the range 
< E{B — V^) < 3 mag, the model was let free to find a fit without imposing any constraint on the star 
formation rate of the various generations. The "best fit" model obtained in this way (CEl) is shown in 
the top panel of Fig. 4 and in Tables 2 and 3: while the shape of the continuum can be fitted with an 
E(B-V)=0.56, this extinction cannot reproduce the equivalent widths of the emission lines and the model 
Ha line is too weak compared to the one observed. 

For comparison, we also present a star burst model with a constant extinction (model CE2, lower 
panel of Fig. 4, Tables 2 and 3): having imposed a current burst|^, "the best fit" model can reproduce 
the continuum with the same extinction as model CEl, but presents far too strong emission lines. In 
conclusion, while a fit to the continuum can be found with an E(B-V) at most ^ 0.6 (as shown by these 
models with constant extinction) , the line ratios necessitate a color excess greater than 1 , as proved by the 
observed Balmer decrement. Selective extinction is required to reproduce simultaneously the strength of 
the lines and the shape of the continuum of e(a) spectra. 



^"This was done imposing that the SFR during the last 10^ was at least twice as much the maximum SFR allowed 
in the old population. 



Table 1: Populations of the starburst models 
Age Model A Model B Model C Model D Model E 





Rv 


= 3.1 


Rv 


= 3.1 


Rv 


= 5 


Rv 


= 3.1 


Rv 


= 3.1 




SF 


E(B-V) 


SF 


E(B-V) 


SF 


E(B-V) 


SF 


E(B-V) 


SF 


E{B-V) 




5.9 


2.60 


91.5 


3.00 


17.2 


2.60 


10.0 


3.40 


0.7 


1.37 


3 • 10^ 


4.7 


1.20 


100.0 


2.41 


18.5 


1.10 


6.8 


2.02 


0.7 


0.93 


8 ■ 10^ 


4.7 


1.03 


95.3 


0.74 


17.8 


0.92 


7.8 


0.72 


0.7 


0.95 


107 


6.4 


1.40 


91.5 


1.57 


15.3 


0.83 


7.5 


1.77 


0.7 


1.04 


5 ■ 10^ 


6.5 


1.17 


93.4 


1.04 


17.3 


1.20 


1.1 


0.70 


0.7 


0.48 


108 


7.2 


0.45 


96.8 


0.72 


12.3 


0.45 


1.0 


0.31 


0.7 


0.22 


3 • 108 


1.2 


0.40 


1.0 


0.02 


3.5 


0.40 


0.3 


0.16 


0.7 


0.23 


5-108 


0.6 


0.40 


1.5 


0.00 


1.7 


0.40 


0.9 


0.50 


0.7 


0.46 


109 


0.6 


0.40 


1.3 


0.03 


1.7 


0.40 


0.4 


0.03 


0.7 


0.24 


old 


1.0 


0.40 


1.0 


0.48 


1.0 


0.20 


1.0 


0.51 


1.0 


0.50 
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Table 2: Populations of the post-starburst (PSB) and constant extinction (CE) models 



Age 


Model PSBl 


Model PSB2 


Model CEl 


Model CE2 


yr 


Ry 


= 3.1 


Ry 


= 3.1 


Ry 


= 3.1 


Ry 


= 3.1 




SF 


E(B-V) 


SF 


E(B-V) 


SF 


E(B-V) 


SF 


E(B-V) 




0.6 


1.02 


0.6 


0.79 


0.05 


0.56 


8.3 


0.56 


3 • 10'^ 


0.6 


0.81 


0.1 


0.46 


0.2 


0.56 


8.3 


0.56 


8- 10® 


0.2 


1.18 


0.1 


0.10 


4.5 


0.56 


8.6 


0.56 


107 


0.1 


0.32 


0.8 


0.70 


4.4 


0.56 


8.3 


0.56 


5 • lO'' 


5.7 


0.65 


5.8 


0.68 


13.6 


0.56 


13.1 


0.56 


10* 


4.7 


0.72 


5.8 


0.68 


1.4 


0.56 


0.1 


0.56 


3 ■ 10* 


2.6 


0.21 


1.6 


0.16 


11.5 


0.56 


22.2 


0.56 


5 ■ 10* 


2.5 


1.22 


1.7 


0.64 


3.8 


0.56 


22.3 


0.56 


109 


2.3 


0.59 


2.1 


0.99 


0.1 


0.56 


18.9 


0.56 


old 


1.0 


0.47 


1.0 


0.52 


1.0 


0.56 


1.0 


0.56 



Table 3: Results 





W(Ha+NII) " 


W(H/3) 


W(H5) 


W([Oii]) 




SFR= 


description 




-62.48±5.3 


0.68±0.8 


5.64±0.5 


-12.68±2.1 


88.0±13.0 




observed e(a) spectrum 








STARBURST MODELS 






Model A 


-65.05 


0.72 


5.52 


-11.40 


28.9 


38 


2 ■ 10* yr burst 


Model B 


-61.26 


0.30 


5.65 


-11.69 


85.0 


102 


2 • 10* yr burst, high FIR/5500 


Model C 


-56.97 


0.72 


5.13 


-12.71 


93.7 


117 


different extinction law 


Model D 


-59.46 


0.43 


5.33 


-12.23 


62.8 


126 


lO'^ yr burst 


Model E 


-59.03 


0.55 


5.72 


-12.24 


8.3 


9 


~ constant SFR 








POST-STARBURST MODELS 






Model PSBl 


-46.50 


0.71 


5.81 


-12.36 


10.1 


4 


post-SB 1 


Model PSB2 


-56.68 


-2.06 


5.25 


-13.05 


9.5 


4 


post-SB 2 






MODELS WITH CONSTANT EXTINCTION 




Model CEl 


-37.64 


-0.44 


4.94 


-15.49 


9.7 


0.2 


Constant ext., free SF 


Model CE2 


-379.54 


-64.47 


-6.97 


-152.86 


12.3 


13 


Constant ext., SB model 



a) All equivalent widths are in A, negative in emission and positive in absorption. The model W(Hq+NII) was found 
adopting Nil = 0.2Ha which is the mean value observed by Tresse et al. (1999) for this Ha strength. 

b) LpjB/Lv=Lpjjj/{5500 x Lssoo) where Lpm (ergs^^) is the total FIR luminosity and Lbboo (ergs~^ A ) 
is the luminosity at 5500 A as derived from the spectrum. 

e) Current star formation rate in M^yr^^ derived from the model fitting the observed luminosity at 5500 A assuming 
a Salpeter IMF (x=2.35, 0.1-100 Mq). 

d) Errors quoted are standard deviations of the bootstrap median value within the sample of 19 galaxies. 
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MODELS WITH CONSTANT EXTINCTION 

20 r 




FIR/V=9.7 



15 r FIR/V=12.3 

" I 1 

: 

10 r 



6 7 8 9 
log(AGE) (yr) 



10 



Fig. 4. Results of the models with constant extinction listed in Table 2. Symbols as in Fig. 2. 



5. Conclusions 



Spectra with a strong B.S line in absorption and a moderate [Oil] line in emission - named e(a) 
spectra - are numerous among luminous infrared galaxies. We have shown that the optical e(a) spectrum 
can be well reproduced by a spectrophotometric model that includes both the stellar component and the 
emission lines and continuum of the ionized gas, assuming dust extinction to vary with the age of the stellar 
population. 

The main conclusions of this work can be summarized as follows: 

1) On the basis of our models, two conditions are necessary to explain the e(a) spectra: 

a) the youngest stellar generations (age < 10^) in e(a) galaxies must be highly obscured by dust. Their 
extinction values are significantly higher than those applying to the previous stellar populations. Models 
with a constant value of extinction are unsuccessful in fitting the optical spectrum; 

b) these galaxies must host strong ongoing star formation at a level higher than in quiescent spirals. 
Post-starburst models with a small residual star formation fail to explain simultaneously the strength of the 
Ha and of the H/3 lines. 

2) The degeneracy of the optical spectrum hinders a quantitative estimate of the burst duration and 
intensity. The spectrum is consistent with, but does not require, a "long" starburst that began 10^ — 10^ yr 
ago. 

3) A model fitting the e(a) spectrum does not necessarily produce the FIR/optical ratio observed: the 
optical spectrum (including the observed Ha luminosity and Hq/H/3 ratio) does not constrain the current 
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star formation rate. The observed FIR/V can be explained by starburst models if a significant fraction of 
the FIR luminosity originates in regions that are practically obscured at optical wavelengths. 

4) All the results listed above refer to the case of a dust screen placed in front of each stellar generation. 
Models in which dust and stars are uniformly mixed fail to fit the e(a) spectrum, producing too low a 
reddening of the emerging emission lines. The observed reddening seems to be at least partly due to 
foreground dust. 

5) The extinction correction factor at Ha found from the Balmer decrement in the e(a) Very Luminous 
Infrared galaxies of our sample {E{}ia) ~ 15, corresponding to a median E(B-V)=1.1, W98a) is in 
agreement with that derived with the same method in other FIR luminous samples (Kim et a. 1998, Veilleux 
et al. 1995, median E(B-V) for non-AGN galaxies =1.13 and 1.05, respectively). This factor is about 6 
times higher than the average E(Ra) = 2.5 in nearby spirals (Kennicutt 1992). Ignoring the fact that the 
slit might have missed a fraction of the Ha luminosity and applying the extinction correction to derive the 
total current star formation rate from the relation SFR(M0yr~^) = 0.8 x 10"''^ E(Ha)) x L(Ha) erg s""'^ 
(Kennicutt 1992)|^, yields a value of the SFR that is still underestimated by a factor ~ 3 compared to the 
SFR derived from the FIR luminosity. 
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